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We perform molecular dynamics simulation of a small number of particles in a box with periodic 
boundary conditions from a view point of chaotic dynamical systems. There is a transition at a 
critical energy Ec that each particle is confined in each unit cell for E < Ec, and the chaotic diffusion 
occurs for E > Ec- We find an anomalous behavior of the jump frequency above the critical energy 
in a two-particle system, which is related with the infinitely alternating stability change of the 
straight motion passing through a saddle point. We find simultaneous jump motions just above the 
^ , critical energy in a four-particle system and sixteen-particle system, which is also related with the 

motion passing through the saddle point. 
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(-H . I. INTRODUCTION 
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^ ■ Liquid-solid phase transitions have been long studied by the molecular dynamics (MD) simulation, since Alder 
^ ' and Wainwright found the liquid-solid phase transition in a system of hard spheres Various thermodynamic 
^jj . properties have been studied with numerical simulations of a large number of particles 0. Alder and Wainwright 
^ ' also found the long-time tail of the velocity auto-correlation, which leads the divergence of the transport coefficient 
c/2 , such as the diffusion constant and the viscosity in one and two dimensions Q. Recently, it is considered that 
dynamical heterogeneities play an important role in the glass transition and the jamming transition of granular 
material 4[. Successive jump motions and intermittent swirling motions are observed in MD simulations of supercooled 
liquids [1, A deterministic chaotic dynamics is assumed in the MD simulation, and therefore, the liquid-solid 

phase transition might be interpreted as a kind of dynamical transition in chaotic systems. The liquid -solid phase 
transition was studied using the Lyapunov spectrum which characterizes the chaotic dynamics by Posch et al. The 
Q most-positive Lyapunov exponent exhibits a maximum at the phase transition Isl 9||. On the other hand, Awazu 
O . and Munakata-Hu studied a system of only two hard disks in a rectangular box [lOl . They showed that there is a 
dynamical transition similar to a liquid-solid phase transition in the system of a small number of particles. . 
I \ In this paper, we study numerically a similar dynamical transition in Hamiltonian systems of a small number of 
^ • particles in a box with periodic boundary conditions. That is, each particle is confined in each unit cell for E < E^, 
■ and the chaotic diffusion occurs for E > Ec- We call it a delocalization transition in this paper. The delocalization 
C transition in a small number of particles has a similarity to the liquid-solid phase transition in an infinite-size system, 

' in that each particle is confined in each unit cell in the solid phase and the particles move around beyond the unit 
, cells in the liquid phase. However, the relation between the delocalization transition and the statistical-mechanical 
T-H ' liquid-solid phase transition is not clear now. The purpose of this paper is to show some peculiar behaviors found 
, near the delocalization transition. 

In §2, we study a two-particle system, and discuss an anomalous behavior of the chaotic diffusion above the critical 
energy. In §3, we study a four-particle system and a sixteen-particle system, and discuss simultaneous jump motions 
just above the critical energy. We will relate these behaviors with the motion passing through a saddle point. 

We consider Hamiltonian systems interacting with repulsive forces. The repulsive force between two particles at ri 
and r2 is assumed to be F(r) = —{dXJ/dr) ■ r/|r| where r = r2 — ri, using the Lennard-Jones potential U{r): 
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for r < Tq with a = ro/(2^/^). The attractive part of the Lennard-Jones potential is not used, that is, F{r) is assumed 
to be for r > tq. The repulsive force is therefore a short-range force, and U{r) is a monotonically decreasing function 
of r. For our numerical simulation, tq is set to be 1, e = 0.01 and the mass m is assumed to be 1 for all particles. 
The fourth order symplectic method of timestep At — 0.0001 was used for the numerical simulation. 
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FIG. 1: Trajectories of the first particle at (a) E = 1.992 and (b) E = 2.228 for L = 1.2. (c) Trajectory of the first particle at 
E — 2.104. The radius of a circle surrounded by the four shaded regions is denoted by R and the gate width is denoted by d. 
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FIG. 2: (a) Potential energy along a trajectory (-L/4, y) of the first particle for L = 1.2. (b) Critical energy of the delocalization 
transition as a function of L. 



II. ANOMALOUS BEHAVIOR OF CHAOTIC DIFFUSION IN A TWO-PARTICLE SYSTEM 

As one of the simplest system, we consider a two-particle system in a square box. The two particles are further 
assumed to be located at point-symmetric positions with respect to the center {L/2, L/2) of the square box. That is, 
the position r2 of the second particle is represented as r2 = (x2, 2/2) = {L — xi, L — yi) by the position ri = (xi ,yi) of 
the first particle. The periodic boundary conditions are further assumed. This system is equivalent to a system that 
the pair of particles are arranged periodically in space with period L. The equations of motion for the first particle 
are given by 
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where {x[,y'j) = {x2 + i ■ L,y2 + j ■ L) and 
particles at {x[,y'j) satisfying rij < vq 
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(2/1 ~ y'j)^- The summation is taken only for the 



This system is a Hamiltonian system with four degrees of freedom of motion 
and conserves the total energy E. In this paper, the potential energy is counted as AC/ = U — U{rQ) from the 
minimum value U{ro). Our system is one of the minimum system that can exhibit the delocalization transition. 
Figures 1(a) and (b) display the trajectories of {xi,yi) for (a) E = 1.992 and (b) E = 2.228 at L = 1.2. The 
initial conditions are slightly different, i.e., (a) xi(0) = yi(0) — L/4:,dxi{0)/dt = 0.05,dyi{0)/dt = 1.4 and (b) 
a;i(0) = 2/1(0) — L/A,dxi{0)/dt — 0.05, dyi(0)/dt — 1.5. Note that the ranges of the x and y coordinates are very 
different. For E = 1.992, the particle exhibits chaotic motion but it is confined in a square box of [0, L/2] x [0, L/2]. 
However, for E — 2.228, the particle exhibits chaotic diffusion. The {x, y) space is divided into cells with unit length 
L/2. The particle is confined inside of one cell when E is smaller than a critical value Ec, however, the particle 
moves around in a chaotic manner from one cell to another cell for E > Ec- This transition is called a delocalization 
transition. 
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FIG. 3: (a) Frequency / that the particle passes through the periodic boundaries as a function of E for L = 1.2. (b) Time 
evolution of yi{t) for two initial velocities dy\{Q)/dt = 1.52 (dashed line) and 1.56 (solid curve), (c) Poincare plot in the xi-Vx 
space at the section yi = L/A. The initial conditions are Vx(Q) ~ 0.005 x n (n — 1,2,3 ■ ■ ■ 10.), Vy{0) — 1.52, x(0) = I//4 and 
y(0) = L/4. 
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FIG. 4: (a) Linear growth rate A of \5xi \ as a function of E for L = 1.2 in the range between 1.75 < E < 3. (b) Linear growth 
rate A in the range between 2.0876 < E < 2.0877. (c) Width AEn of the nth parameter range of E where A is nearly 0. The 
dashed line denotes a line of AEn oc 1/(2.4)". 



The point (I//4, L/4) is a point where the potential energy takes the minimum. On the other hand, the point 
(L/4, L/2) is a saddle point of the potential energy, because the potential U{L/A,y) takes a maximum value at 
y = L/2 along a line x = L/4 but U{x, L/2) takes a minimum value at x = L/4 along a line y — L/2. Figure 2(a) 
displays the potential energy AU{y) = U{ri) — //(rg) along a line (L/4, y). The minimum energy to get over the saddle 
point is therefore calculated as Ec = 2{C7(L/2) — f7(ro)}, because particles located at (3L/4, L/2) and (—L/4, L/2) 
interact with the first particle at the saddle point (L/4, L/2). It is a critical energy for the delocalization transition. 
Figure 2(b) displays the critical energy E^. as a function of the system size L. As L is decreased, the critical energy 
Ec increases rapidly, because the two particles are confined in a smaller box. 

The shaded region in Fig. 1(a) is surrounded by four arcs whose radius R satisfies 2{L''(2i?) — U{ro)} — E. As 
E is increased, R decreases, because J7 is a monotonically decreasing function. For E < Ec, R is larger than L/4. 
Then, the shaded region surrounded by the four arcs is confined in the cell of size L/2, and the chaotic motion is also 
confined inside of the unit cell region. At E = Ec, i? is equal to L/4 and the four arcs touch the saddle points at 




FIG. 5: (a) Period t2 of yi{t) in Eq. (3) as a function of In \E — Ec\. (b) Stability exponent A' as a function of the period t2 for 
a simple model Eq. (5) . 
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y = L /2 or X = L /2 and the chaotic motion invades into the neighboring cells, and the chaotic diffusion appears. For 
E > Ec, the shaded region surrounded by the four arcs is connected with the shaded region in the neighboring cell 
as shown in Fig. 1(c) at E' = 2.104. The arcs of the neighboring four cells form a circle of radius R. The gate width 
d formed by the overlap of the neighboring shaded regions at the section y — L/2 is equal to d — L/2 ~ 2R, which is 
proportional to E — Ec near the critical energy Ec. The particle moves to the neighboring cells through these narrow 
gates with gate width d ior E > Ec- 

We have numerically calculated the frequency / that the first particle passes through the periodic boundaries 
xi = UxL or yi — UyL {rix and Uy are integers). The initial condition is a;i(0) = yi{0) — L / 4, dxi{0) / dt = 0.001 
and the initial velocity dyi{0)/dt is changed as a parameter to change the total energy E. Figure 3(a) displays the 
frequency / as a function of the energy E for L — 1.2. The frequency / is for i? < ii^^ ~ 2.09, because the particle 
is confined inside of the cell as shown in Fig. 1(a). For E > E^ f{E) takes a nonzero value, however, the behavior of 
f{E) is rather anomalous. Roughly speaking, there are two branches for f{E), and which branch is selected depends 
on the energy. In the lower branch of f{E), f{E) increases as / oc \E — Ec\ near the critical energy, which corresponds 
to the chaotic diffusion as shown in Fig. 1(b). It is because the frequency / is proportional to the gate width d formed 
by the overlap of the neighboring shaded regions, and therefore / is proportional to E — Ec near the critical energy 
Ec. 

In the energy intervals where f{E) takes the upper branches in Fig. 3(a), the particle exhibits an almost straight 
motion in the y direction. This mode of motion is called the accelerator mode. This type accelerator modes were 
studied in various systems such as the standard map or the chaotic diffusion in an oscillatory Benard convection ^12..|13l. 
The anomalous chaotic diffusion in a dissipative system was also studied in an oscillatory Benard convection [14 1. 
Figure 3(b) displays time evolutions of yi{t) for two initial velocities dyi{0)/dt = 1.52 and 1.56. The other initial 
values are the same as before, i.e., a;i(0) — yi{0) = L/4,dxi{0)/dt = 0.001. For Vy{0) = 1.52, yi{t) increases in 
proportion to t, which implies the straight motion in the y-direction. For Vy{0) = 1.56, yi{t) exhibits chaotic motion 
around j/i = 0. Figure 3(c) displays the Poincare plot in the x-v^ space at the section yi = L/4 for the initial 
conditions Vx{0) = 0.005 x n (n is an integer), Vy{0) = 1.52, x(0) = L/4 and y(0) = L/4. A straight motion of 
xi{t) = L/4 =const. corresponds to a central fixed point in this Poincare map. The KAM tori around this fixed 
point represent stable accelerator modes. A large sea of chaos exists around the small region of stable KAM tori. 
If the initial conditions are in the KAM region, the particle exhibits a straight motion, and if the initial conditions 
are outside of the KAM region, the particle exhibits chaotic diffusion. If the initial conditions are located in the 
intermediate region, the straight motion appears intermittently and the diffusion constant becomes large as shown in 
[12], [13] and [14]. 

The completely straight motion along the line x = L/4 is a special solution. The straight motion is stable for some 
parameter ranges where the KAM tori appear as shown in Fig. 3(c). The stability of the straight motion can be 
studied from the time evolution of the small deviation Sx(t) from the straight motion along the line x = L/4. The 
y-coordinate yi{t) of the straight motion and the deviation 5x{t) obey 




(3) 



dt^ 
d'^Sxi 

I,] K ^ - - - ^ ■ J -J 

(4) 

The motion of yi{t) is periodic in time and the period increases to infinity at the critical energy Ec. The transverse 
instability is measured by the linear growth rate A of |i5a;i(i:)|. Figure 4(a) displays A as a function of E in the range 
1.75 < E < 3. The unstable parameter regions and the stable regions appear alternatively. The unstable parameter 
regions for E > Ec correspond to the parameter regions where the lower branch of f{E) appears in Fig. 3(a) and the 
stable parameter regions for E > Ec correspond to the parameter regions of the upper branch of f{E) in Fig. 3(a). 
The width AE' of each parameter region decreases as E approaches Ec. Figure 4(b) displays the magnification of 
Fig. 4(a) in the range 2.0876 < E < 2.0877. The alternation of the stability continues infinitely as E ^ Ec both 
for E < Ec and E > Ec- Figure 4(c) displays the width A_E„ of the nth stable parameter ranges as a function of n 
in a semi-logarithmic scale. The width AE„ for E < Ec is marked by + and the one for E > Ec is marked by x , 
although they are almost overlapped in Fig. 4(c). Figure 4(c) implies that the width A_E„ decreases exponentially as 
AE (X (1/2.4)" as E approaches Ec. 

To the best of our knowledge, this type of singular behavior was not reported before. This behavior can be 



5 



qualitatively understood with a simpler model equation. The equation for dxi is roughly approximated as 

(fSxi 



+ LofSxi = 0, for < t < ti, 
+ uilSxi = 0, for ti<t<t2, (5) 



where t2 is the period of yi by Eq. (3), t2 — ti is the time during which the particle stays near the saddle point 
yi = L/2, and ti is the time during which the particle is far away from the saddle point. That is, the complex 
function of the left-hand side of the equation for Sxi in Eq. (4) is approximated at a piocowise linear equation (5), 
which characterizes respectively the dynamical behavior near and far from the saddle point. As E is close to Ec, the 
period t2 increases in a logarithmic manner as t2 ~ —aln\E — Ec\- The approximation by the second equation in 
Eq. (5) would be better near the critical energy Ec, because t2 — ti increases infinitely and the particle stays near 
the saddle point for a very long time. Figure 5(a) displays t2 as a function of In |_E — Ec\ by the direct numerical 
simulation of Eq. (3). The parameter a is evaluated as a ~ 0.0827. The parameter a is related to the behavior of 
the potential AU near the saddle point shown in Fig. 2(a), i.e., AU Ec — l/(2a^)(j/ — L/2)^. The stability of the 
linear equation (5) per one period t2 can be calculated from eigenvalues of the matrix C = AB where 



A 



cosa;2^2 sina;2^2 

-(0)2/0)1) Sina;2t2 {U!2 / U)]) cos U!2t2 



Q_( COS wi^i cosa;2ii + (oii/a;2) sinwiii sina;2ti sinwiti cos W2^i — (oii/a;2) cos wiii sinw2ii 
~ \^ coswifi sina;2ii — (t^i/'^2) sinwiti cosa;2ii smujiti sina;2ii + (wi/0'2) coswiti cosuj2ti 

The largest eigenvalue A' of C expresses the stability of the linear equation (5) per one period t2- Figure 5(b) displays 
the largest eigenvalue A' if it is a real number, or the modulus of A' if the eigenvalues arc complex number, as a function 
of t2 for Lol ~ 100 and wl = 1940.7. The growth of the deviation |5a;i(t)| in a large time interval t is evaluated as 
|5a;i(t)| <^ |A'|*^*^ ^ exp(At), because the particle passes through the saddle point {t/t2) times in the time interval t. 
Therefore, the exponent A for the transverse instability is evaluated as A ^ In |A'|/i2- Here, the value = 1940.7 is 
evaluated from the right-hand side of Eq. (4) at the saddle point (i/4, L/2) near E = Ec, and loI is roughly estimated 
from the right-hand side of Eq. (4) near the fixed point (L/4, L/4). If |A'| > 1, 5xi increases to infinity and the 
straight motion is unstable. Figure 5(b) shows that the stability alternates periodically as a function of the period t2- 
The period T of the alternation of the stability is equal to 77/0)2 ~ 0.0713, because of the periodicity of the matrix 
A. (The sign of A' has no meaning.) The period t2 of the straight motion is evaluated as ^2 — — alnl-En — Ec\ at 
the nth energy interval. It is equalized to nT, because it corresponds to the nth stable parameter range of A'. Then, 
\En — Ec\ = exp(— nT/a) ^ {T ja)"'^ ~ (1/2.37)". This exponent 2.37 is consistent with the numerically obtained 
value 2.4 in Fig. 4(c). Note that the exponent is determined only by a and L02, which are further determined only by 
d'^AU/dx^ and d'^AUjdy^ at the saddle point. Thus, we find that the simple approximation by Eq. (5) explains the 
singular behavior rather well. 



III. SIMULTANEOUS JUMP MOTIONS OF A SMALL NUMBER OF PARTICLES 



A. Four particles in a square box 



The anomalous behavior in the two-particle system near the delocalization transition is due to the coexistence 
of the KAM tori and chaos. If the degree of the freedom of motion is increased, the regions of the KAM tori are 
expected to become smaller, and therefore, the effect of the anomalous behavior is also expected to become smaller. 
However, there is another peculiar behavior in chaotic diffusion by a number of particles near the delocalization 
transition. We have performed numerical simulations of four particles in a square box oi L = 1.4 under the periodic 
boundary conditions. No special symmetry is further assumed. Figure 6(a) is a trajectory (x\,y\) of the first particle 
at = 0.89. The chaotic diffusion occurs in the (1, 1) or (1, —1) direction. The unit cell is a square of size L/B in 
this system. Figure 6(b) displays time evolutions of Xi{€) for i = 1,2,3 and 4 for S = 0.89. Jump motions toward 
the neighboring cells are clearly seen. The jump motion occurs almost simultaneously for the four particles. The first 
particle comes into the unit cell of the second particle, then, the second particle needs to go away from the original 
unit cell and moves into the unit cell of the third particle. Successively, the third particle comes into the unit cell of 
the fourth cell particles, and the fourth particle invades into the unit cell of the first particle. We call this type of 
successive jump motions a simultaneous jump motion in this paper. It reminds of the heterogeneous and intermittent 
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FIG. 6: (a) Trajectory of the first particle at _E = 0.89 in a four-particle system in a square box of L = 1.4. (b) Time evolutions 
of x^{t) for i = 1, 2, 3 and 4 at S = 0.89. 
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FIG. 7: (a) Potential energy AU as a function of the simultaneous displacement by r of the four particles along the (1, 1) or 
(1, —1) direction, (b) Displacement of the four particles for the potential AU in Fig. 7(a). Four particles move in the directions 
of the arrows, (c) Potential energy AU as a function of the single-particle displacement by r along the (0, ~1) direction, (d) 
Numerically evaluated diffusion constant as a function of E. 



motion of particles near the glass transition or the jamming transition. The jump motion occurs as a swirling motion, 
because the energy E is close to the critical energy Ec- 

Figure 7(a) displays the potential energy AU, when the four particles move simultaneously as (a;i,yi) = (3L/4-t- 
r,3L/4- r),{x2,y2) - (2L/4 - r, i/4 + r), (xg, yg) = {2L/4 + r,3L/4 + r) and (x4,y4) = {L-r,L/4-r). Here, 
the potential energy AU is calculated as the sum of the potential energy of the four particles measured from U{ro). 
When r = L/16, the potential AU takes the maximum. The point is a saddle point because it is a minimum point 
for the potential energy along the transverse direction. When r = L/8, the potential energy takes the same value as 
the initial configuration at r = 0, because the initial and the final configurations of the four particles are the same 
parallelogram, although their directions are different as shown in Fig. 7(b). Each particle moves by one unit cell 
by making the course shown in Fig. 7(b). The configuration at r = L/16 is considered to be a saddle point. The 
minimum energy to go over the saddle point is i5 = 0.8643. In direct numerical simulation, the chaotic diffusion 
appeared for E > 0.875. The critical energy for the delocalization transition is close to the energy of the saddle point. 
If the total energy E is close to the energy of the saddle point, it is necessary for the four particles to exhibit the 
successive motion as shown in Fig. 7(b) to go over the saddle point. It takes much larger energy for a single particle 
to go over the barrier, setting the other particles fixed in the original position. Figure 7(c) displays the change of 
the potential energy AU for such a course (xi,?/i) = (3L/4,3L/4— r), (x2,y2) = (2^/4, L/4), {x^^ys) ~ {2L/A,3L/4:) 
and {x4,y4) = {L,L/A). When only the first particle is moved downward by L = 1.4, an original configuration is 
recovered owing to the periodic boundary conditions. However, the critical energy is about 1475. Such a jump motion 
cannot occur if the energy E is around 1. Figure 7(d) displays the diffusion constant D for one particle calculated 
from {{x{t) — x{0))'^ + {y{t) — y(0))^)/(4i). The diffusion constant increases roughly in proportion to E — E^. 



B. Sixteen particles in a rectangular box 

Next, we consider sixteen particles in a rectangular box of size x Ly = L x {\/?>/2)L under the periodic boundary 
conditions. The size L is fixed to be 3.008. A triangular lattice is a natural configuration for the solid phase in 
this system. Figure 8(a) displays trajectories of the sixteen particles, which form a typical triangular lattice for 
E = 3.827. Figure 8(b) displays trajectories of the sixteen particles dX E = 8.158 between 800 < t < 1000. String-like 
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FIG. 8: (a) Trajectories of sixteen particles a,t E — 3.827 between 800 < t < 1000. (b) Trajectories of sixteen particles at 
E = 8.158 between 800 < t < 1000. (c) Time evolution of y^{t) for i = 4, 7, 12, 15 and 16 at S = 8.158. 




FIG. 9: (a) Time evolution of pt{l, r) for I = 0.55 and r = 4 at £ = 8.158. (b) q{l, r) as a function of t for Z = L/8 = 0.375 at 
several E^s. (c) X4,il,''') as a function of r for I = 0.375 at several E's. 



trajectories appear near x = 2.5, which imply successive jump motions. Figure 8(c) displays time evolutions of j/i(t) 
for z = 4, 7, 12, 15 and 16. Particles i — 4,7, 12, 16 jump downward by one layer almost simultaneously, and the other 
particles are confined in the original cells in this time interval. The jump motion occurs like a chain-reaction in the 
2/-direction owing to the periodic boundary conditions. This is another example of the simultaneous jump motions. In 
this sixteen-particle system, the simultaneous jump motion by all the sixteen particles rarely occurs near the critical 
energy in contrast to the four-particle system. 

We try to characterize the jump motions using some quantities such as q(l,T) and X4(^;^ which were used to 
characterize the heterogeneous dynamics near the glass transition or the jamming transition [15|, , [l3|. Figure 9(a) 
displays time evolution of pt(Z,T) = J^^iii^ii'^) ^ ^li^ ^ '''))^ + ~ Vii'^ ~~ ~ ^\ for ^ = 0-55 and r = 4 at 

E — 8.158, where Q(x) is the Heaviside step function. The quantity pt{l,T) represents the number of particles whose 
displacement in the time interval r is larger than I. Figure 9(a) implies that the simultaneous jump motion by four 
particles occurs intermittently. The four particles are not fixed but are chosen randomly. Simultaneous jump motions 
by eight particles or sixteen particles do not appear in the time sequence. Figure 9(b) displays q(l,T) = {qt{l,T)) as 
a function of r for several values of E ior I = L/8, where qt{l,T) = (1/16) ^ {(^i(0 ~ ~ '''))^ + iVii'^) ~ 

?/i(t — r))^}^/^] and (• • •) implies the time average. The quantity qt{l, t) represents the number ratio of particles whose 
displacement in the time interval t is smaller than I. The quantity q{l, r) is the temporal average of qt{l, t). The time 
scale, where q{l,T) decreases rapidly from 1 to 0, increases as E approaches Ec- The quantity q{l,T) characterizes 
that the chaotic jump motions become more intermittent as E ^ Ec- Figure 9(c) displays Xii^T) a function of t 
for several values of E for / — L/8, where X4:{1,t) is defined as 

X4il,r)^il/16)[{qt{l,Tf)~{qtil,r)n 

The quantity Xii^T) take a maximum value at a certain time, and the time scale increases as E is decreased to Ec- 
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FIG. 10: (a) Potential energy AUi (solid curve) along a zigzag displacement of four particles, At/2 along a horizontal displace- 
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ment. (b) Diffusion constants Dx and Dy as a function of E. (c) and Dy as a function of E. (d) / ' as a function of 

E, where / is the frequency of the jump motion. 

This quantity also characterizes the intermittency of the simultaneous jump motions. 

Figure 10(a) displays the potential energy AUi (solid curve) along the zigzag displacement of four particles: 
(x4,2/4) = {7L/8 + r/2,V3/2{L/8-r)),{x7,y7) - (6L/8 + r/2, V3/2(3L/8-r)), (xi2, 2/12) = (7L/8-r/2, v^/2(5i/8- 
r)), (ccie, yie) — (i — r-/2.V3/2(7i/8 — r)), and the potential energy AJ72 (dashed curve) along the horizontal displace- 
ment (xi, yi) = (L/8 - r, (V3/2)L/8), (2:2, y2) = (3L/8 - r, (V3/2)L/8), (xg, ya) = (5L/8 - r, (V3/2)i/8), (x4, 2/4) - 
(7i/8 - r, (V3/2)i/8). The other particles are fixed to the original positions. The zigzag displacement corresponds 
the simultaneous jump motion of four particles shown in Fig. 8(b). The potential energy at the peak position is smaller 
for AC/2 than for AUi. It implies that the horizontal motion is easier to occur. However, the zigzag displacement 
appears more often in direct numerical simulation as shown in Fig. 8(b). We do not understand the reason well 
yet, however, it is partly because the configuration number of the zigzag displacement is larger than the horizontal 
displacement. That is, the configuration number of the horizontal displacements is 4 because one of the four layers 
moves, however, the combination number of the zigzag displacement is counted as 4 x 2 x 2 = 16. 

Figure 10(b) displays two diffusion constants, i.e., in the x-direction and Dy in the y-direction, as a function of 
E. The diffusion constants Dx and Dy are slightly different in this numerical simulation because of the anisotropy 

in this system, however, their behavior is very similar near the critical energy. Figure 10(c) displays Dx^^ and D^^^ 
as a function of E. The diffusion constants increase roughly as D^.y ^ [E — Eq)^ for E > 6.8 where Eq ~ 6. When 
E < 6.8, the jump motion becomes very rare, however the jump frequency is not zero for E > Ec ^ 5.2. Here, Ec is 
estimated by numerical simulation performed until t = 100000, and it might be lowered in a even longer simulation. 
The critical energy for the delocalization transition is comparable to the maximum energy of AU shown in Fig. 10(a). 
Figure 10(d) displays f{Ey/^, where f{E) is the frequency of the jump motion per unit time, as a function of E. The 
frequency / of the jump motion also increases as / ~ (i5 — Eq)^ for E > 6.8, although the reason of the exponent 3 
is not understood yet. 

IV. SUMMARY AND DISCUSSION 

We have performed numerical simulations of Hamiltonian systems with repulsive interaction, and studied the 
delocalization transition. We have found that the dynamics around the saddle point is very important near the 
delocalization transition, because the saddle point plays a role of the lowest barrier between the neighboring cells. In 
a two-particle system, we have found infinite alternation of the stability and the instability of the straight motion 
passing through the saddle point. It is closely related to the anomalous behavior of the jump frequency / above the 
critical energy. In a four-particle system in a square box and a sixteen-particle system in a rectangular box, we have 
found that simultaneous jump motions occur intermittently near the transition point. If each particle is located deeply 
inside of each cell, the chaotic motion for the particle is weakly correlated from that for the other particles. However, 
it is necessary for several particles to exhibit jump motion almost simultaneously to go over the saddle point using 
a small amount of energy. All the four particles exhibit the simultaneous jump motion in the four-particle system. 
On the other hand, four particles exhibit the simultaneous jump motion near the delocalization transition in the 
sixteen-particle system. These simulatenous jump motions might be related with the heterogeneous dynamics near 
the glass transition and the jamming transition, although the detailed relation with such systems of a large number 
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of particles is not known well. 

There is a definite critical energy of the delocalization transition in Hamiltonian systems of a small number of 
particles, below which each particle is completely confined in each unit cell. The critical energy seems to be determined 
by the energy of the lowest saddle point. We consider that the critical energy of the delocalization transition in the 
Hamiltonian systems of a small number of particles is generally different from the critical energy for the statistical- 
mechanical liquid-solid phase transition in an infinitely large system. It is partly because a long-range order needs 
to appear against the statistical-mechanical fluctuations below the liquid-solid phase transition point. For example, 
the long-range order does not appear in two dimensions owing to the long-wavelength fluctuations [18J. The long- 
wavelength fluctuations are not directly related to the energy of the saddle points. However, the detailed relation of 
the delocalization transition and the statistical-mechanical liquid-solid phase transition is not understood yet but it 
is left as a future problem. 
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